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ABSTRACT OF THE DISSERTATION 

Observation of High Frequency Internal Waves 
in the Pacific Ocean 

by 

Robert Libman Zalkan 

Doctor of Philosophy in Oceanography 
University of California, San Diego, 1968 
Professor Fred N. Spiess, Chairman 



Observations of high frequency internal waves were 
made from 1 to 14 December 1966, in the deep sea off 3a Ja 
California. The fluctuations of the depth of an isotherm 
were measured with a three-element horizontal array 
attached to a stable platform. 

The waves are characterized as a broad band 
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phenomenon with a continuous distribution in frequency. For 
short Intervals of time, a narrow frequency band within the 
continuum is adequately described as a horizontally plane 
wave of a single vertical mode. Furthermore, this simpli- 
fied structure is stationary over time spans of several 
days. High modes are present in the low frequency waves. 
Above 4 cycles per hour, however, the first mode predomi- 
nates. 

The spectral shape is consistent with the shear 
limited equilibrium spectrum proposed by Phillips (1966). 

In addition, the predominance of the first mode at high 
frequencies further emphasizes the importance of shear 
instability in Internal wave propagation. 

The horizontal properties of the wave field indi- 
cate well-defined directions of narrow-beamed propagation. 
These directions and the dispersive properties of the 
propagation have led to the identification of local topo- 
graphic features as generating areas of internal waves. 
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CHAPTER I 



INTRODUCTION 

Since internal waves constantly change the interior 
structure and velocity fields of the ocean, they may have 
important effects on dynamic height calculations, current 
measurement, mixing, marine organisms, underwater sound pro- 
pagation, and submarine operations. An adequate description 
of the wave field is required in order to determine the ex- 
tent of these effects. In addition, an understanding of the 
mechanisms of .generation and decay is necessary to the 
determination of the role internal waves play in energy 
transfer processes. Specifically, the description should 
Include modal structure, energy spectra, and directional 
information. 

At the present level of theoretical development, a 
model that adequately describes these features is clearly 
not feasible, and little is known about the mechanisms of 
generation and decay. Our purpose, therefore, is to seek 
such a description from observation of high frequency 
internal waves in the deep sea. 

The complex structure of Internal waves makes both 
data collection and analysis difficult. Unlike surface 
waves, internal waves do not have a unique dispersion 
relationship. Instead, there is a different relationship 
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for each member of the Infinite set of possible modes. In 
addition, the maximum vertical displacements of these waves 
occur at intermediate depths rather than at the free sur- 
face, and are, therefore, difficult to observe. 

To overcome these difficulties, simultaneous measure- 
ments of the fluctuating depth of an Isotherm were recorded 
with a three-element horizontal array. This procedure 
allows the establishment of a phase difference (travel time) 
relationship leading to a vector wave number and, therefore, 
directional determination. The wave number magnitudes can 
then be tested against the theoretical dispersion curves for 
confirmation of progressive wave phenomena. 

Stokes (1847) first showed that internal waves can 
exist in a stratified fluid consisting of two homogeneous 
layers. Later Webb (1884) and Greenhill (1887) dealt with 
the multi-layered case. In 1933, FJeldstad extended the 
theory to Include a continuous density gradient. More 
recent systematic treatments of the subject are by Eckart 
(I960) and Tolstoy (1963). A summary of the theory of small 
amplitude internal waves in still water is given in 
Appendix A. 

Several investigators have attempted to obtain di- 
rectional information by towing a thermistor chain in a 
sequence of traverses of different direction. This pro- 
cedure allows one to obtain a directional description of 
the wave field by Interpretation of measured Doppler shifts. 
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Using this method, Metcalf et al (1962) reported internal 
waves of 17 minute periods traveling south from the equator. 
However, LaFond (1964) in the Pacific and Charnock (1965) 
in the Atlantic, had little success. 

Ufford (1947) appears to be the first to make 
simultaneous measurements of internal waves at spatially 
separated points. His first method was to make repeated 
lowerings of bathythermographs from the bow, midsection, 
and stern, of a ship located over the San Diego trough, in 
1200 meters of water. The results were inconclusive, due 
to the small phase differences and, therefore, large error 
limits in the determination of wave number. To increase the 
phase differences, Ufford moved to shallow water (40 meters) 
and employed both bathythermograph lowerings from three 
anchored ships and a buoyed array of three thermocouples. 
Unfortunately, each observation spanned only a few hours of 
time, and he used only those portions of the records 
(average 1 hour, 7 minutes) in which a single frequency 
predominated. 

Further work in shallow water (18 meters) has 
been accomplished at the Navy Electronics Laboratory’s 
Oceanographic Tower. Here, a three-element array of iso- 
therm followers has been successfully used to determine 
directional information (LaFond 1962a, Cox 1962). Neverthe- 
less, mid-oceanic observations are still required since 
shallow water propagation differs considerably from that in 



deep water. In particular, the wave form, and therefore the 
frequency content, Is affected by the nearness of the sea 
floor (LaFond, 1961b) and the propagation direction by shore- 
ward refraction (LaFond, 1962b ) . Also, there is probably a 
preference for the first mode at all frequencies, since the 
shallow depth, and thus narrow waveguide, means large shear 
for modes greater than one. 

We, therefore, seek our information from the deep 
ocean. The difficulty of employing the array technique in 
deep water is attested to by the previous absence of such 
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CHAPTER II 



OBSERVATIONAL EQUIPMENT 

A three element horizontal array was obtained by 
attaching heavy duty booms to the Research Platform 
FLIP's (Fisher & Spiess, 1963) bow or upper structure. The 
attachment was done by hinges enabling the booms to be 
properly stowed while FLIP was in the horizontal position. 
Cables to the sea-sensing units were fed over the ends of 
the booms. Upon "flipping" to the vertical position, the 
booms were rigged out to form an equilateral triangle of 
30.5 meters (100 feet) per side (Figure 1). Each boom had 
its own Independent sensing system consisting of 1) an 
Isotherm follower, 2) a vertical staff of four variable 
frequency thermometers, and 3) a vibrotron pressure guage. 

The isotherm follower is controlled by a thermistor 
in a sensing unit. A servo system is arranged to vary the 
depth of the sensing unit in such a way that the temperature 
measured by the thermistor remains constant. The thermistor 
signal is developed across a DC bridge. One leg of the 
bridge is connected by a cable to a thermistor in the sea- 
sensing element. The other leg is composed of a variable 
bank of resistors for selection of the desired temperature. 
When the thermistor has the same resistance as the variable 
resistors, the bridge is in balance. When the temperature 
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Figure 1. Triangular array attached to the Research Platform FLIP. 



at the thermistor changes, the thermistor changes resistance 
and a DC voltage occurs across the bridge. This voltage Is 
fed through magnetic amplifiers and silicon controlled 
rectifiers to the windings of a DC motor. The polarity and 
amplitude of the signal determine the direction and speed at 
which the motor turns. The motor drives a winch which 
raises or lowers the sea-sensing element until the bridge 
Is again In balance. In this way the thermistor In the sea- 
sensing unit Is kept at the depth of an Isotherm and hence 
the name Isotherm Follower (LaFond, 1961a; Good, 1968). 

The upper part of the sea-sensing element Is a 
hollow pressure case. A vlbrotron pressure guage was Instal- 
led In this portion and exposed to sea pressure via an oil- 
filled chamber with a flexible diaphragm. The output of the 
vlbrotron was used for determining the depth of the sea- 
sensing unit. 

Since servo-systems are subject to dead bands (i.e., 
a certain level of signal Is necessary before the system 
responds), a staff with four variable frequency thermometers 
(Squler, 1967) was Installed on each sea-sensing unit to be 
used as a vernier. The thermometers were placed .6 meters 
(2 feet) and 1.8 meters (6 feet) above and below the servo- 
system thermistor. A scheme of interpolation was then used 
to determine the actual depth of the isotherm (see Appendix 
B). 

The data were recorded both digitally for later 
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analysis and by analog methods for on the spot monitoring. 

The digital data collecting system consisted of 1) a clock, 

2) a sequential sampling unit, 3) a counting system, and 4) 
an incremental tape recorder. Every twenty seconds the 
clock put out a pulse which triggered the sampling unit. 

The sampling unit consisted of a stepping switch which se- 
quentially connected each of the channels of information to 
the counting system for one second's duration. The counting 
system was a period counter set to determine the number of 
microseconds required for 4096 cycles of signal. This 
number, along with a serial time count, was then recorded on 
magnetic tape by an incremental recorder. The information 
recorded from each of the sea-sensing units consisted of 
the vibrotron pressure signal, the four variable frequency 
thermometer temperatures, and the vibrotron pressure signal 
again. Each follower was sampled in turn and then the 
azimuthal orientation of the array was recorded to conclude 
the cycle. After the 19 pieces of information had been 
recorded, there was one second of dead time. Then the 
clock sent out a pulse and the whole procedure started again. 

DC voltages for analog recording were obtained by 
sending the three vibrotron signals and one temperature 
signal to frequency meters. The output voltage from the 
meters drove the pen displacement system of a four channel 
chart recorder. 

Since FLIP is free to rotate about the vertical 
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axis, It was necessary to employ an azimuthal orientation 
system. An error signal between the desired heading and 
FLIP's actual heading and a rate of rotation signal from the 
gyrocompass repeater were used to drive the tilt boxes of 
two hydraulic motors. The motors in turn drove two off-axis 
propellers. The torque from the propellers would rotate 
FLIP until the desired heading was achieved. If FLIP 
drifted off the desired heading the system would again be 
activated to maintain the desired heading. 

A block diagram of the data collection and orient- 
ation systems is shown in Figure 2. 

Preliminary data were obtained during February, 
March, and April of 1966 while Flip was operating in 
Hawaiian waters. Analysis of these data indicated the need 
to : 

1. increase the sensitivity and reliability of the iso- 
therm followers 

2. incorporate the thermometer staffs as verniers 

3* change the data recording system to incorporate 2 

4. reposition the booms and winches 

5. change the cable fairleads because of 4 

6. replace lightweight telescoping booms with heavy-duty 
ones 

7* employ the orientation system to keep the array stabil- 
ized in azimuth for as long a time as possible. 

The implementation of these modifications resulted 
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Figure 2. Block diagram of data collection and orientation systems. 
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In the system described above. A short seatrlal was con- 
ducted near Catalina Island In November. Then FLIP was 
towed to deep water off Baja California and data were 
collected during December 1966. 



CHAPTER III 



OBSERVATIONS 

A station southwest of San Diego (near 30°N,121°W) 
was maintained from 1 December 1966 to 14 December 1966. 

FLIP was allowed to drift with the mean Integrated current 
over her 300 foot draft. Since the waves travel relative 
to the current, this drift minimized any Doppler effects. 

A track of FLIP'S drift is shown In Figure 3. 

Due to the aerodynamic characteristics of FLIP's 
upper structure, she tends to orient with her keel Into the 
wind. Once she has self-oriented, the propeller system has 
little trouble maintaining the heading. If the wind 
direction shifts more than twenty degrees, however, the load 
becomes too great on the propeller system and re-orientation 
is required. It was possible, nevertheless, to break the 
data into 15 short series during which the variance of 
array orientation was small. The duration, mean azimuth of 
the keel, and azimuthal standard deviation are given in 
Table I. 

Time records of series 4, 5, and 6 are presented in 
Figure 4, and are typical examples of the internal oscil- 
lations during the entire sampling period. Each series 
contains isotherm depth vs. time from all three followers. 
The middle trace corresponds to the depth scale on the left. 
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Figure 3. Chart of FLIP’S track and location 

of generating areas (after H. W. Menard, 
I960), Numbers on track refer to the 
day in December 1966 (0800 position). 
Letters refer to generating areas. 
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TABLE I 

Division of Data into 15 Series 



•TIME (Dec. 1966) 


MEAN HEADING (°T) 


STANDARD 
DEVIATION ( 


061400-021120 


329.2 


4.8 


021120-031005 


239.8 


3.2 


021827-031712 


240.2 


3.1 


031816-041701 


240.7 


2.6 


041701-051547 


240.6 


1.2 


051819-061704 


239.1 


2.8 


070115-071633 


348.5 


2.7 


071633-080625 


328.5 


4.1 


080751-090636 


343.2 


4.3 


090636-100515 


342.4 


1.9 


100557-110443 


011.8 


2.3 


101430-111315 


010.4 


3.5 


120338-122108 


002,4 


27.3 


122108-132003 


356.5 


6.8 


132003-140556 


007.6 


10.6 
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Fisure 4. Depth of the observed Isotherm vs time. 

The upper and lower traces have been offset + 10 meters. 



The upper and lower traces have been offset + 10 meters 
respectively. The interpolation scheme used to calculate 
the isotherm depth may be found in Appendix B. 

The first characteristic one should notice in the 
time records is the predominance of the tidal oscillations. 
Tidal frequencies contain most of the energy and are, there- 
fore, responsible for the largest vertical motions of the 
Isotherm. 

Superimposed on the tidal motions are the high fre- 
quency fluctuations with which this thesis is concerned. 
Several Interesting properties of the high frequency 
motions are apparent in the time records. First, the 
records from the three elements exhibit good visual coher- 
ence throughout the data. Second, the larger amplitudes 
of the high frequency oscillations seem to come and go as 
if they were modulated. And finally, the oscillations 
appear to be always present, indicating, at least in a 
qualitative manner, that they may be stationary. These 
conclusions are Justified quantitatively in Chapter V. 

The second property, i.e. that of amplitude modu- 
lation, has been observed in shallow water by LaFond (1961) 
and may be caused by two effects. The most likely cause of 
this amplitude fluctuation is the beating effect. When 
energy is present as a continuous spectrum, beating occurs 
among the various frequencies and a rapidly varying ampli- 
tude envelope is developed. 
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The second possible cause is a resonance condition 
described by Eckart (1961) and hereafter referred to as 
Eckart Resonance. When the waves are propagated in a well- 
localized group — localized in space, frequency, and wave- 
number — a resonant coupling between the layers of high 
stability can occur. The effect of this resonance is to 
cause the energy to oscillate between the layers. If the 
conditions are variable along the propagation path, and one 
measures the motion in one layer only, an amplitude modu- 
lation will be observed, with the modulating frequency 
dependent upon the sharpness of the layers and their verti- 
cal separation. The same effect would be observed if the 
source is not stationary in a spatial sense. That the 
necessary conditions for this resonance existed is shown in 
Chapter IV. 



CHAPTER IV 



ANALYSIS TECHNIQUES 



The power spectrum has become a basic tool for pre- 
senting the properties of long complicated time series in a 
condensed form. Here, the method of first Fourier trans- 
forming and then summing over the squares has been employed. 

Consider a discretely-sampled, real series of finite 
length, X(J), J = 0, N-l. The X(j) may be represented by 
the complex series 



X(j)=Z*(m)e l27rim/N - j= o , N-l 

m = o 



(l) 



( 2 ) 



Then the complex Fourier coefficients are: 

X(m)=J E X(j) e i27rim/N , m = 0 , N-l 

j-° 

Taking into account the symmetry of the transform of a real 
series, a spectral estimate of elementary frequency bands 
with two degrees of freedom is obtained by 



S x (m)= A m X(m) X(m) , m = 0 , ^ - 1 



(3) 



where * denotes complex conjunction and Aq = 1, m=o 

Am. = 2, m>o 

To increase spectral stability, a sum over adjacent ele- 
mentary bands is formed. 
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q 

» Ajl a 

S x (h c ) = Z S x (m) = £ A m X (m) X(m) 

m=p 



(4) 



where p and q are the lower and upper band Indices over 

which It Is desired to sum and h = p + q . Letting 

c ©n u ©r 2 

a = q - p + 1 (the number of bands Included In the sum), 

equation ( 4 ) yields 2a degrees of freedom per estimate . 

Since It Is desirable to model a continuous speotrum of 

Infinite extent which has dimension length 2 /cycle/time = 

2 

L T, It Is necessary to divide S(h) by the band width, a. 



Thus 

E x (h c ) = J^Am X*(m) X(m) 



and simple summing is replaced by band averaging. 

The conversion from the Indices m, h c , a, to fre- 
quency Is accomplished by using the sampling Interval A t. 



Let NAt = T 

then lim =2ir(jAt)f 

N hQ 

so that f=~ » ^center - J~ » an dAf — "[* 



A further discussion of these techniques Is pre- 
sented by Tukey In Harris (19 66). Haubrlch (1965) shows 
that these methods produce the same result as those of 
transforming the correlation function. The actual calcu- 
lations were performed on a CDC 3600 digital computer using 



the fast Fourier transform techniques of Cooley and Tukey 
(1965). 



Once the series has been transformed, the cross 
spectral estimate between two series is simply 




(6) 



The estimates of coherence and phase difference 
are respectively: 

|G xy (h c ) 



R xy (h c ) = 



E x (h c )E y (h c ) 



(7) 



and 

<£ xy (h c ) = Arg (G xy (h c )) y y leads x 

Considering the triangular array 

3 I 

I ^;(h c )=Z Arg(G,(h c )), j. 12,23, 31 

i a 




(8) 



(9) 



For a single train of plane waves, equation (7) 13 
identically equal to 1 , and equation (9) is identically 
equal to 0. Loss of coherence and a non-zero may be 

caused by a finite beam width (Munk et al, 1963; Cox, 1962). 



Cox also discusses reduction of coherence caused by turbu- 
lence and variable phase velocity. It should be noted that 
If a = 1 (l.e. p = q) equation (7) has only one degree of 



to zero. Therefore It Is necessary to band average. 

Failure to average over at least two bands would lead to 

the conclusion of plane wave propagation, when some other 

form of propagation may exist. 

Using the three phase differences defined In 

equation (8), It Is possible to determine three vector wave 

numbers. Let vector wave number k have components (k , k ). 

x y 

The three phase differences, for plane wave propagation 



tween elements of the array. Since there are three 
equations and only two unknowns, (k x , ky), the wave number 
Is over determined. Considering the equations two at a 
time, three separate solutions for the vector components 
may be found. 




become 




( 10 ) 



where X^> and Yj are the spatial separation components be- 



Then 




( 11 ) 




where 



m = 1,3 



(12) 
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i and J being unit vectors in the x and y directions and m 
indexing the three sets of solutions. 

The mean direction of propagation is defined as 



8 = Arg (k, + k 2 + k 3 ) 



(13) 



Beam width indication,^© , is defined as 



A 9 = max(Arg(ki))-min(Arg(kj) i,j=l,3;i*j a*) 



where the branch cut is taken to form^© as the minimum 
angle between vector solutions k^ and kj (as in figure 5). 




Figure 5 

Definition of a e 

It should be noted that if Rj = 1 and, thus,£(J^ = 0, 
then only two of the equations are independent, the k m 's 
are Identical, and^6 = 0. These results agree with the 
earlier discussion of plane wave propagation. 



23 



Once a vector wave number with a small A 9 has been 
found, it becomes necessary to determine whether the motion 
is due to internal waves or to some other mechanism such as 
anisotropic turbulence. The distinction between wave motion 
and motions due to other causes can be made with the help of 
diagnostic dispersion curves. The solutions , k rma are, 
plotted on an (J~k diagram and tested against the theoreti- 
cal dispersion curves to see if there is agreement. It must 
be remembered, however, that solutions which fail the test 
may still be wave motions not due to single plane waves (i.e. 
multiple modes, or two or more plane waves of the same mode 
and frequency traveling in different directions, etc.). 

To obtain the dispersion curves, it is first neces- 
sary to calculate density, P (z), from the observed 
temperature, T(z), and salinity, S(z), distributions. Next 
p(z) is used to calculate N(z) . Finally, for each fre- 
quency band of interest, the wave equation is numerically 
integrated to solve for the k n . 

Figures 6 and 7 show the observed temperature and 
salinity distributions, the calculated density and stability 
frequency distributions, and the theoretical dispersion 
curves for the first 10 modes. Figure 6 reflects conditions 
Just prior to the commencement of wave data collection, and 
Figure 7 reflects the conditions immediately after termi- 
nating data collection. The temperature distribution in the 
upper 200 m was obtained by manually lowering an isotherm 
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TEMPERATURE ( C) 

S 10 15 20 

SALINITY 10/001 STABILITY FREQUENCY (CPH) 

33.20 33. TO 3H.20 0 6 12 18 






Figure 6. Propagation conditions for 1 December 1966 
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Figure 7. Propagation conditions for 14 December 1966 
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follower. The distributions of salinity and of temperature 
below 200 m were obtained from nearby CalCOFI stations 
(Figure 6 in December, Figure 7 in May). The difference in 
the two salinity distributions at shallow depths is probably 
attributable to sampling at different seasons as opposed to 
being a permanent spatial feature. In any case, temperature 
is the controlling factor within the thermocline and the 
salinity distributions were used for "order of magnitude" 
purposes. Since spatial aliasing will occur for X ' > .016, 
the dispersion curves are terminated at this value. Below 
400 meters T(z) and N(z) monotonically decrease while S(z) 
and p (z) monotonically increase with increasing depth. The 
values of these parameters corresponding to figures 6 and 7 
are tabulated in Table II. 

Several features of Figures 6 and 7 are worth 
further consideration. Both stability frequency distri- 
butions have many relative maxima. These maxima are 
reflected in the dispersion curves as wiggles (Fig. 6, 3rd. 
mode) or changes in curvature (Fig 6, 2nd mode). It is the 
coupling between these layers of high stability which pro- 
duces Eckart Resonance. The greater the vertical 
separation between the layers, the sharper the resonance. 

The resonance regions are shown in the dispersion curves as 
areas where two modal lines converge and then diverge. For 
example, in Figure 7, a resonance region between the second 
and third mode is apparent at a frequency slightly greater 
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than 4cph (X -l == .0085)* Th© coupling seems to occur be- 
tween the layers centered about 75 meters and 250 meters. 
Other regions of resonance are also present. 

TABLE II 



Distribution of Temperature, Salinity, Density, 
and Stability Frequency below 400. 

1 December 1966 



Depth Temperature 

500 5.93 
1000 3.91 
1500 2.76 
2000 2.10 
3000 1.63 
4000 1.56 



500 5.92 
1000 3.96 
1500 2.89 
2000 2.12 
3000 1.65 
4000 1 . 58 



Salinity 


Density 


34.20 


1.02993 


34.47 


1.03320 


34.57 


1.03622 


34.63 


1.03912 


34.66 


1.04472 


34.66 


1.05021 



14 December 1966 



34.22 


1.02990 


34.49 


1.03321 


34.58 


1.03622 


34.64 


1.03913 


34.68 


1 . 04474 


34.69 


1.05024 



Stability Frequ ency 



2.28 

2.06 

1.93 

1.77 

1.76 



3.00 

2.28 

2.06 

1.93 

1.80 

1.77 



The two sets of stability frequency and dispersion 
curves show several Interesting differences. The upper 
high stability layer Is sharper and slightly deeper In 
Figure 7 relative to figure 6. On the other hand, the 
lower layer is weaker and much deeper. Increasing the 
separation between layers and causing the 2nd - 3rd mode 
resonant region to shift to a lower frequency and inverse 
wavenumber. While the first modal curve shifted only 
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slightly and only in the low frequency region, all of the 
higher modal curves shifted toward lower frequencies. The 
weakening of the Intermediate relative maxima may also ex- 
plain why the curves of Figure 7 are much smoother than 
those of Figure 6. Notice, however, that the regions of 
slightly negative curvature still remain in modes 2 and 3. 

These differences may have been caused by both 
temporal and spatial changes in the internal structure. 

But, because the internal waves are continually changing 
the subsurface structure, it is possible that the mean con- 
ditions were stationary throughout the data collection 
interval# and that the differences between figures 6 and 
7 were caused by taking the two profiles when different 
phase relationships existed among the waves. In this case 
these figures would represent perturbations about the mean 
conditions at the two spatially and temporally separated 
points. 

Figure 8 shows the mean depth of the 16.85°C iso- 
therm for series 1-12 and the 16.00°C Isotherm for series 
13-15 at each of the three sea-sensing elements. The 
16.85°C isotherm shows large oscillations during the first 
12 series. For the last three series, however, it was 
necessary to shift to a lower temperature in order to re-- 
maln within the thermocllne. Therefore, it is probable that 
the differences are due to a combination of both mechanisms. 
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Figure 8. Mean depth of the observed Isotherm vs time. 



CHAPTER V 



RESULTS OF ANALYSIS 

The fifteen aeries tabulated in Table I were 
analyzed by the methods of Chapter IV. As examples, the 
results for series 1, 4, 8, 10, and 14 are shown graphically 
in Figures 10 through 14. The layout of each figure is as 
follows: 

The upper left hand graph shows the three energy 
spectra - one for each isotherm follower. Spectral averag- 
ing (equation (5) ) was done over eight elementary 
frequency bands yielding 16 degrees of freedom per spectral 

estimate. For 16 degrees of freedom, there is 80$ confi- 

2 2 2 2 2 
dence that .58cr <S <1.47cr and 95$ confidence that .43cr <-S 

2 P 2 

<1.80<r where S is the spectral estimate and c f is the 

true spectral value. Continuing down the left hand side of 
the page, the second graph shows the three pair coherences 
as defined in equation (7). For 16 degrees of freedom, the 
95$ significance level for coherence is .25. The third 
graphs display ^ <J> as defined in equation (9). The upper 
right hand graph shows the computed wave number magnitudes 
(equation (11) ) plotted as points and the theoretical dis- 
persion curves (modes 1 to 10) of 1 December 1966 plotted 
as lines. Since the modal lines converge in the lower fre- 
quency region, an expanded scale plot for series 1 is shown 



30 



INVERSE WAVELENGTH (1 /METERS) 

.-0005 .0010 .0015 .0020 .0025 .0030 .0035 



31 




FREQUENCY (CPH) 



Figure 9. Expanded scale dispersion diagram for series 1. 



SUM PH, ANG. (RAO I RMS) COHERENCE (R**2) LOG ENERGY PENS ITT (M**2/CPH) 
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FIGURE 10 SERIES I 

ENERGY SPECTRA? PAIR COHERENCES? 
AND 

SINGLE WAVE SOLUTIONS 
011400-021120 DECEMBER 1966 



SUM PM. ANG. (RAO IANS) COHERENCE <R**2) LOG ENERGY DENSITY <M»*2/CPH) 
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FIGURE II SERIES 4 



ENERGY SPECTRA? PAIR COHERENCES? 
AND 

SINGLE WAVE SOLUTIONS 
031816-041701 DECEMBER 1966 



SUM PH. BNG. (RROIRNS) COHERENCE (R**2) LOG ENERGY DENSITY (M*»2/CPH) 
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FIGURE 12 SERIES 8 



ENERGY SPECTRA? PAIR COHERENCES? 
AND 

SINGLE WAVE SOLUTIONS 
071633-080625 DECEMBER 1966 



SUM PH. ANG. (RAO IANS) COHERENCE (R**2) LOG ENERGY DENSITY j-2/CPHI 
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FIGURE 13 SERIES 10 

ENERGY SPECTRA? PAIR COHERENCES? 
AND 

SINGLE WAVE SOLUTIONS 
090636-100515 DECEMBER 1966 



SUM PH. ANG. (RADIANS) COHERENCE (R**2) LOG ENERGY DENSITY (M**2/CPH) 
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FIGURE 14 SERIES 14 

ENERGY SPECTRA? PAIR COHERENCES? 
AND 

SINGLE WAVE SOLUTIONS 
122108-131954 DECEMBER 1966 
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in Figure 9 (corresponds to boxed area of figure 10) . The 
final graph displays the directional information. Each line 
runs from (©--^-^) to (©+— -^) where © and A© are defined by 
equations (13) and (14) respectively. 

Various narrow band spectral peaks are apparent 
throughout the 15 series. However, the peak heights lie 
between the 80$ and 95$ confidence limits and the fre- 
quencies at which they occur change from series to series. 
The peaks are caused by nearby frequencies having the 
proper phase relationships to accentuate their particular 
frequency band during the short duration of sampling. This 
result is consistent with the expected results from sampling 
a continuous, random, Gaussian process. 

The coherence is quite high for the lower frequen- 
cies but in some cases falls off quite rapidly in the 5 to 
7 cph region. Similarly, £ <f> remains close to zero when the 
coherence is high and then appears to oscillate randomly 
as the coherence drops. The rate at which the coherence 
falls is greatest when the energy spectrum is also falling 
rapidly with frequency and vice versa. A possible expla- 
nation for the fall off is that in these cases there 
appears to be a lack of first mode solutions at the higher 
frequencies. For modes greater than 1, the wavelength 
becomes smaller than the array spacing near these frequen- 
cies. Therefore as frequency increases more and more 
wavelengths occur between the elements of the array (see 
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Cox 1962). In other cases, however, the coherence rises 
again to quite high values and then falls off in the 10 cph 
region with £ $ "behaving in a consistent manner. In these 
cases, there appears to "be a predominance of first mode 
solutions above 5 cph (compare series 1 and 8). 

Some of these solutions, particularly above 5 cph, 
lie to the right of the first mode dispersion curve but 
have good coherence, small X <j>, and narrow beam width, A0. 
In several series, these solutions actually seem to form a 
separate modal line below that of the first mode. The best 
example can be found in series 8 where there is a multitude 
of narrow beamed high frequency waves traveling toward 270°. 
The most likely cause of this phenomenon is a Doppler shift 
due to current shear. FLIP drifts with the mean integrated 
current over her entire 300 foot draft and these high fre- 
quency waves are superimposed on the mean integrated current 
over their relatively narrow high stability layer (see Fig- 
ures 6 and 7). The difference between these currents 
required to account for the observed Doppler shifts is on 
the order of only 5 cm/sec. The details of the necessary 
calculations and a specific example may be found in 
Appendix C. 

The solutions considered acceptable as single wave 
propagation, together with their modal numbers and mean 
propagation directions are listed in Table III. Solutions 
falling below mode 1 are included as mode 1 and those 
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TABLE III 



Summary of Acceptable Solutions 



frequency 


Band No. 
CPH 


2 

.33 


3 

.68 


4 


5 


6 


7 


0 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


32 


33 


34 


35 


36 


37 

12.63 




1.03 


1.38 


1.74 


2.09 


2.44 


2.79 


3.14 


3.49 


3.84 


4.20 


4.55 


4.90 


5.25 


5.60 


5.96 


6.31 


6.66 


7.01 


7.36 


7.71 


8 . 06 


8.42 


8.77 


9.12 


9.47 


9.82 


10.17 


10.52 


10.88 


11.23 


11.58 


11.93 


12.28 


SERIES 


Mode 


5 


10 


4 


2 


c; 


3 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 






1 










1 


1 




1 




















X 


Direction 


101 


082 


033 


035 


077 


102 


153 


. 118 


006 


055 


348 


093 


086 


112 


102 


126 






108 










303 


268 




332 




















o 


Mode 


10 


4 


4 


2 


8 


2 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 
































d 


Direction 


194 


281 


312 


021 


165 


047 


152 


128 


028 


111 


010 


083 


115 


347 


136 


075 


178 


138 


150 


128 


157 
































—7 


Mode 


A 


3 


3 


1 


5 


1 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


1 
















3 


Direction 


340 


307 


233 


039 


07 5 


085 


154 


148 


224 


033 


078 


052 


089 


131 


138 


132 


192 


161 


127 


136 


165 


231 


215 


168 


254 


221 




235 


220 
















i 


Mode 


10 


9 


5 


3 


2 


2 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 




1 


1 


1 
















4 


Direction 


233 


170 


132 


036 


130 


105 


010 


045 


137 


059 


023 


055 


006 


020 


240 


275 


286 


184 


145 


249 


341 


354 


185 




318 




184 


280 


339 


















Mode 


10 


4 


5 


2 


2 


3 


2 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


1 


1 




















5 


Direction 


141 


210 


095 


094 


110 


114 


352 


122 


069 


062 


295 


121 


309 


334 


315 


185 


170 


166 


247 


251 


276 


263 


290 




285 


286 


278 






















Mode 


5 


3 


3 


9 


3 


3 


1 


2 


1 


1 


1 


1 


1 




l 


1 


1 


1 


1 


1 


1 


1 


1 




1 








1 
















6 


Direction 


225 


046 


194 


151 


233 


187 


060 


105 


154 


115 


107 


121 


330 




021 


276 


263 


345 


269 


257 


293 


282 


046 




272 








295 
















7 


Mode 

Direction 


10 

262 


10 

322 


10 

348 


4 

306 


10 

194 


2 

.112 


2 

323 


2 

232 


1 

262 


3 

285 


2 

019 


1 

012 


2 

215 


1 

047 


1 

210 


1 

315 


1 

250 


1 

270 


2 

240 


2 

032 


2 

059 


1 

252 


1 

299 


1 

276 


1 

255 


1 

26l 


1 

264 


1 

256 


1 

279 


1 

276 


1 

315 


1 

292 










8 


Mode 

Direction 


10 

082 


10 

301 


9 

325 


10 

224 


6 

027 


2 

198 


3 

290 


3 

344 


3 

201 


2 

177 


1 

276 


1 

249 


2 

261 


2 

318 


1 

121 


1 

221 




1 

233 


1 

231 


1 

241 


1 

313 


1 

269 


1 

276 


1 

268 


1 

2 50 


1 

276 


1 

27 1 


1 

260 


1 

272 


1 

277 


1 

273 


1 

293 


1 

269 




1 

268 


1 

272 


9 


Mode 

Direction 


9 

317 


10 

146 


10 

296 


6 

276 


10 

265 


4 

245 


4 

275 


2 

271 


1 

117 


1 

194 


1 

113 


2 

214 


1 

343 


1 

248 


1 

225 




1 

005 


1 

274 


1 

265 


1 

265 


1 

267 


1 

255 


1 

280 


1 

001 


1 

275 


1 

241 


1 

271 






1 

257 


1 

271 












10 


Mode 

Direction 


10 

266 


10 

260 


10 

327 


10 

291 


10 

358 


3 

002 


5 

345 


3 

051 


2 

286 


2 

327 


1 

321 


2 

346 


1 

336 


1 

331 


1 

062 


1 

056 


1 

227 




1 

255 






1 

310 


1 

282 


1 

141 


1 

113 










1 

010 






1 

094 








j 11 


Mode 

Direction 


10 

059 


10 

340 


10 

334 


10 

336 


3 

192 


2 

342 


3 

038 


2 

307 


3 

355 


2 

328 




2 

316 


2 

010 


1 

057 




1 

032 


1 

053 




1 

054 


1 

020 


1 

001 






1 

307 


1 

308 


1 

320 


1 

338 


1 

054 


1 

037 




1 

324 


1 

170 




1 

018 


1 

271 


1 

097 


12 


Mode 

Direction 


10 

270 


10 

357 


10 

358 


4 

299 


2 

277 


4 

049 


2 

350 


3 

270 


3 

Oil 


2 

322 


1 

003 


1 

278 


2 

006 


1 

054 


1 

024 


1 

108 




1 

031 


1 

065 


1 

065 


1 

000 


1 

147 






1 

310 






1 

007 




1 

274 


1 

339 




1 

239 








13 


Mode 

Direction 


10 

269 


6 

025 


7 

075 


5 

035 


4 

349 


3 

038 


4 

025 


1 

358 


2 

041 


1 

013 


1 

298 


1 

036 


1 

030 


1 

008 


1 

017 


1 

006 




1 

041 


1 

030 




1 

053 


1 

002 




1 

002 


1 

007 


1 

008 








1 

005 














14 


Mode 

Direction 


10 

036 


6 

286 


9 

180 


7 

032 


3 

109 


1 

015 


1 

051 


1 

042 


1 

120 


2 

315 


1 

166 


1 

204 


1 

319 


1 

220 


11 

190 


1 

245 


1 

273 


1 

285 


1 

023 


1 

042 






1 

278 






















1 

077 






15 


Mode 

Direction 


10 

301 


10 

269 


10 

032 


1 

354 


5 

354 


2 

206 


2 

047 


3 

269 


6 

224 


2 

004 


1 

322 


1 

260 


1 

279 


1 

004 


2 

035 


1 

198 


1 

349 


1 

245 


1 

087 


1 

034 


1 

031 


1 

055 


1 

060 


1 

186 


1 

353 




1 

300 




1 

020 




















40 



falling above mode 10 are included as mode 10. All other 
modal indexes refer to the nearest modal line. The 
solutions were also tested against the 14 December dis- 
persion curves. Complete modal agreement was found except 
at 4.9 cph in series 8. Here, the 1 December curves indi- 
cate third mode and 14 December curves indicate second mode 
but at a 2nd-3rd mode resonant region. 

Although high modes frequently occur at the low 
frequencies, the first mode predominates in the high fre- 
quency range. This feature is probably caused by the wave 
guide narrowing with increasing frequency - thus subjecting 
the higher modes to greater shear until they eventually 
become unstable. In general, it is felt that the modal 
solutions are quite good. Notice the good agreement between 
the wave number solutions and the theoretical modal lines in 
Figure 9. Credence to this belief is also given by the good 
fit obtained in the slowness plots of the next chapter. It 
must be remembered however that any particular modal so- 
lution, and especially those with large wave numbers, may 
have been shifted relative to the dispersion curves by the 
Doppler effect discussed above. Table IV lists the number 
of solutions in each mode. 
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TABLE IV 



Number of Acceptable Solutions In Each Mode 



Mode 


Number of Solutions 


1 


257 


2 


47 


3 


23 


4 


11 


5 


9 


6 


5 


7 


2 


8 


2 


9 


5 


10 




Total 


392 



Percent of Solutions 
66 
12 
6 
3 
2 
1 

1/2 

1/2 

1 

8 

100 # 



Numerous wave solutions which appear to be station- 
ary over several days may be found by referring to Table III. 
For example the solutions for f c = 2.44 cph show a 2nd mode 
Internal wave traveling 153 - 1° during series 1, 2, and 3. 
Similarly, for f Q = 7.71, the solutions for series 5, 6, 7, 

8, and 9 show © = 265 - 15° and for t Q = 8.06, 6= 279 ± 3° 
during series 8, 9, and 10. 

To define the general shape of the energy spectrum, 
the 45 spectra (three for each series) were ensemble 
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averaged. Considering that the three spectra from 
each series are not independent, and allowing for the over- 
lapping in time of several of the series, the resulting 
spectral estimates have 194 degrees of freedom. For 194 

degrees of freedom there is 80$ confidence that .87 o 2 <S 2 

2 2 2 2 
*1.138 and 95$ confidence that .8lcr *1.21cf . A log-log 

plot of the average spectrum is shown in Figure 15. 

The spectrum has been broken into characteristic 
regions according to slope. The slope of the curve and the 
power of frequency at which the spectrum falls are identical. 
Initially (.2 cph to 1.7 cph) where the wave guide ( OJ < [\J (^ ) ) 
occupies the entire water column, the spectrum falls as 

-2 2 i 

f ’ . At 1.7 cph, the slope changes to f ** x and the wave 
guide begins to narrow with succeeding spectral estimates 
(i.e., for = 1.74 cph, A 2 = 4000 m for f c = 2.09 cph, 
Az = 1430 m , and for f c = 2.44 cph. A 2 =860 m ; where A z is 
the thickness of the wave guide). From 5 cph to 12 cph, 
the spectrum falls as f "5.25^ Above 12 cph, no internal 
waves appear to be present (the 4 acceptable solutions 
above 12 cph in Table III have been Doppler shifted ), and 
the rate of fall of the spectrum begins to decrease until 
at 50 cph it reaches the inherent noise level of measure- 
ment. This region of zero slope indicates a root mean 
square noise amplitude, £ rms (h c ), of .014 meters. 

The height of the peak at 4.2 cph lies between the 
80$ and 95$ confidence limits. Also, the individual series 
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FREQUENCY (CPH) 



Figure 15. Average energy spectrum. 



solutions for this frequency do not indicate any apparent 
reason for the peak's existence. Therefore, it is felt 
that the peak has no statistical significance. 

The region between 1.7 and 5 cph coincides with the 
region of good wave solutions and the energy appears to be 
attributable to Internal waves. Using a Richardson's 
number Instability argument, Phillips (1966) has proposed a 
theoretical equilibrium spectrum for first mode propagation 
within the thermocllne for frequencies such that the stabil- 
ity layer is small compared to the total depth. His 
spectrum is proportional to f for kd « 1 and to f 
for kd » 1 (where d is the thermocllne depth). Since pre- 
ceding analysis has shown many solutions of mode higher than 
one within the region kd « 1, the agreement between the ob- 
served spectrum and the equilibrium spectrum Indicates that 
the f power law may apply in general for kd<< 1, and is 
not limited solely to the first mode. The change in slope 
from f to f "3.25 occurs where there is a rapid narrow- 
ing of the wave guide, and, as Phillips predicted, where 
kd = 1 (see Figures 6 and 7). For kd » 1, the first mode 
is predominant and the agreement between the observed 
spectrum ( f "3*25) an d the equilibrium spectrum ( f 
furnishes additional support to Phillips' Theory. 



CHAPTER VI 



DIRECTIONAL PROPAGATION AND POSSIBLE SOURCES 



To determine whether there were any preferred 

directions of propagation, the accepted solutions (Table 

III) were subjected to further analysis. Each solution was 

spread out over its beam width and normalized so that 
- 360 ° 

J S(^) d v// = I (15) 

where S(*/ / ) = I/A 9 j 9 -A 9/2 - — 9 + A 9/2 

SW = 0 i t<6-A9/2, V>9+A9/2 

The solutions were then quantized in one degree intervals 
and added together. 



P(j)= z S:(j) ; j = 1/2, 3/2, ... 719/2 
i 



( 16 ) 



where i indexes the solutions and J represents the mid- 
direction of the interval. Finally a 10 degree wide run- 
ning mean filter was used to smooth the results. 

i+f 

QU)=I P(mU=l,2...360 ; P(-n) = P(360-n) <w) 

m=f-§ 



The resulting beam patterns are shown in Figure 16. The 
top pattern Includes all fifteen series. The remaining 5 
patterns Include three series each and show the changes 
which occur with time. Several strong lobes are noticeable 
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Figure 16. Horizontal beam patterns 



4 ? 



within the patterns. The prominent lobe toward 270° pro- 
vides an interesting example. It is virtually absent during 
the first 3 series. It begins to appear during the next 3 
series, and for series 7 to 9, overshadows any other 
direction. It then decreases but remains present for the 
final 6 series. The existence of these lobes leads to the 
conclusion that preferential directions indeed do exist. 

The existence of directional propagation leads, in 
turn, to the consideration of possible sources. Azimuthal 
information is provided by the beam patterns. To determine 
range, the dispersive character of the wave propagation is 
utilized. The wave energy travels at the group velocity, 



which is the slope of the local tangent to the dispersion 



distance of travel is the same for all wave groups, then for 
two wave groups arriving at t and t + ^ t: 




(where n indexes the mode) 



curves at ( f J/X n ) .If the source is localized, so that the 



X = G(f 0 ) t 
X = G(f,)(t+At) 



( 18 ) 

(19) 



so that 



XG' 1 (f 0 )= 1 
XG"' (f,)=t+At 



( 20 ) 

( 21 ) 



Letting AG _, = G _l (f,)-G _, (fo) and subtracting equation (20) 
from equation (21) X = At/AG" 1 



( 22 ) 
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Thus, rang© to a source is obtained by making slowness plots 
(time of arrival vs. inverse group velocity). The slope and 
t-intercept of a reasonable straight line fit give range and 
time of generation respectively. The slowness plots for the 
lobes centered at 270°, 070° and 025° are shown in Figures 
17, 18 and 19. The slope for 270° (= 1.59 day-meters/sec 

gives a distance of approximately 75 miles. Looking in the 
reciprocal direction (l.e., 090°) at a distance of 75 miles 
from FLIP's track, one observes a prominent ridge of sea- 
mounts (labeled "A" in Figure 3). 

In Figure 18 several different slopes are apparent. 
The steepest slope (7.14 day-meters/sec) indicates a distance 
of 330 miles and again, the source seems to be a cluster of 
seamounts (labeled "B" in Figure 3). The next steepest 
slope (2.18 day-meters/sec) indicates a distance of 100 
miles and the seamount marked M C" in Figure 3 seems to be 
the source. The least slope (.92 day-meter/sec) indicates 
a distance of 43 miles. Since this solution occurs during 
the time when FLIP was rapidly drifting toward the south, it 
is felt that the seamount at "D" (Figure 3) was the source 
and that the wave group also was subjected to advection in a 
southerly direction. 

In Figure 19, only the slope of the heavy line 
(5.00 day-meters/sec) corresponds to a topographic feature 
(range = 233 miles - area labeled "E" in Figure 4). The 
other lines give ranges of 110, 345, and 24 miles. They are 
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INVERSE GROUP VELOCITY (sec/meter) 

Figure 17. Slowness plot for direction 270°. Numbers refer 
to frequency band numbers of Table III. Super- 
scripts indicate mode. 



TIME (DAYS -December, 1966) 
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Figure 18. Slowness plot for direction 070°. Numbers refer 
to frequency band numbers of Table III. Super- 
scripts indicate mode* 



TIME (DAYS -December, 1966) 
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INVERSE GROUP VELOCITY (sec/meter) 



Figure 19. Slowness plot for direction 025°. Numbers refer 
to frequency band numbers of Table III. Super- 
scripts Indicate mode. 



shown to Indicate that many of the straight line fits give 
ranges to which no correlation with topography was found. 

The quantity of solutions used to determine the 
slopes makes the range determinations subject to doubt. 

These results were presented not as absolute proof that 
seamounts act as Internal wave sources, but merely to 
demonstrate that the solutions point to them as likely pos- 
sibilities. This conclusion is supported by Sandstrom 
(1966) who has shown that topographic features tend to 
increase the energy In internal waves. 



CHAPTER VII 



CONCLUSIONS 

High frequency internal waves in the deep ocean are 
characterized as a broad band phenomenon with a continuous 
distribution in frequency. The spectral shape is consistent 
with Phillips' theory of a shear limiting process, and the 
predominance of the first mode at high frequency gives 
additional emphasis to the importance of shear instability 
in Internal wave propagation. 

Over short intervals of time, a narrow frequency 
band within the continuum displays great simplicity in 
modal and directional structure and can be adequately de- 
scribed as a horizontally plane wave (l.e., small beam 
width) of a single vertical mode. In addition, it has been 
shown that this simplified structure is statistically 
stationary for time-spans on the order of several days. 

This structural simplicity indicates that the waves are well 
dispersed and Implies the existence of fairly localized 
(both in time and space) sources. 

Considering the properties of the total (broad band) 
wave field, additional evidence of localized sources is 
found. The horizontal nature of the wave field displays 
definite anisotropy, indicating well-defined directions of 
narrow-beamed wave propagation. While the exact character 
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of the anisotropy may change appreciably In times of the 
order of several days, well defined directions of narrow- 
beamed propagation always appear to be present. 

The time-history of the wave field displays the 
presence of dispersive wave trains. As well as Indicating 
localized sources, their presence makes possible the calcu- 
lation of source range. These range determinations, coupled 
with the well-defined propagation directions, Indicate that 
nearby topographic features such as seamounts act as gener- 
ators of high frequency Internal waves. 

In view of the wide distribution of seamounts 
throughout the oceanic basin, the superposition of arrivals 
from the multitude of sources should produce a confused 
"internal sea". Therefore, the waves must not persist over 



long propagation paths 
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APPENDIX A 



Theory of Small Amplitude Internal Waves 

I. List of Symbols 

C velocity of sound 

f Inertial frequency 

9 acceleration of gravity 

k wave number magnitude 

k* x-component of wave number vector 
ky y-component of wave number vector 
n mode Index 

t time 

X horizontal coordinate 

y horizontal coordinate orthogonal to x 

Z vertical coordinate (positive up) 

D depth of bottom 

F force 

N stability frequency 

q specific volume 

£ vertical displacement 

p density 

(jj angular frequency 

dV Infinitesimal volume 

AZ small vertical displacement 

/•*\ 

( ) second derivative with respect to time 
( )-£2 second derivative with respect to Z 
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II. Infinitesimal Waves 

The following presentation consists of only those 
aspects of Internal wave theory necessary to the understand- 
ing of the observational and analytical techniques employed 
In this dissertation. More detailed formulations may be 
found In Eckart (I960) and Tolstoy (1963). 

Consider a vertically stratified, compressible fluid 
at rest In a gravity field. If a parcel of fluid Is raised 
vertically from its equilibrium position and then released, 
it will tend to sink with the force: 

F=g(yo-(/D+^ Az) dV)-yo-3- 2 AzdV (ia) 

the first term being due to buoyancy and the second compres- 
sibility* The vertical equation of motion becomes: 

/ »dVz--g|£zdV-/»§IzdV (2A , 



For infinitesimal motions, the solution to equation (2) is 
that of a harmonic oscillator with the frequency: 




Equation (3A) defines the Brunt-Vaisala or stability fre- 
quency. In the upper ooean, the compressibility term is 
generally small oompared to the buoyancy term and is usually 
neglected. Thus we have harmonic motion caused by buoyancy 
forces due to density stratification. 



Starting with the three-dimensional equations of 



motion, mass conservation, and invoking the following 
assumptions : 

1) N = N(z) and is independent of horizontal 
coordinates 

2) Mean velocities are equal to zero 

3) Boussenesq Approximation - Density variation 
may be neglected except where it causes buoyant forces 

4) Nonlinear terms are negligible 

we arrive at the first order equation of motion 



N 2 (z)-w 2 

c zz -^7 ~ v h£ = 0 



( 4A) 



where OJ is the angular frequency and f is the inertial fre- 
quency due to rotation of the earth and £ is the vertical 
displacement from equilibrium. 

If we assume a solution of plane horizontally pro- 



gressive waves of the form 

,i(k x x + k y y-u>t) 



C(x,y,z,t) = C(z)6 



(5A) 



and k 2 “ k 2 + k 2 



then (4A) becomes 



£ 



zz 



+ 



N 2 -ou 2 

(x) z - f 2 



k 2 £ = 0 



( 6A) 



Internal waves have their maximum amplitude below 
the surface of the water column and at this maximum 



6 0 



C zz <o 



but 



kz =-k ; 

£ 



N 



CU 



2 -o , 2 

2_ f 2 



Therefore either f ^OU orN^cU^f in the oceanN^^>f^ 

so that the former condition Is the correct one, and f and 

N max are the lowest and highest frequencies permitted to 

free Internal waves of form (5A). 

2 £ 2 

If CU 25> ] then equation (6A) becomes 

N 2 -^ 2 

+ k C = 0 



(7A) 



ZZ 



The boundary conditions for our purposes are 

£ = 0 at Z = 0 surface 

£ = 0 at z=-D bottom 

The free surface boundary condition may be used at the sur- 
face but has very little effect on the Internal modes. 

For arbitrary N(z) and fixed CU , solutions to 
equation (7A) are found by Iterative numerical techniques - 
starting with £ = 0 at z =-D and working to satisfy the sur- 
face boundary condition. For a given CU , solutions can be 
found for an Infinite number of discrete values of K. For 
each such value of K^, there Is a corresponding amplitude 
function £ n . The Index n relates to the number of ampli- 
tude maxima within the water column, and (n-1) to the number 
of nodes or zero crossings between the maxima. 

Solutions will be progressive ifCU< Nfz) somewhere in 
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the water column. £ ^ ) 

Cl)<N(z) and "exponential" 



will be "trigonometric" where 
where U)>N (z) (Bee Figure 1A)» 



N C 





Figure 1A. 2nd Mode Amplitude Function 



Within the thermocline where the density gradient, 
and therefore N(z), Is a maximum, the appropriate vari- 
ation of physical parameters and their effects on specific 
volume are given in Table IA 



TABLE IA 



The Variations of Physical Parameters in the 
thermocline and their effects on specific volume 



Parameter 

Temperature 

Salinity 



Variation 
-8° C 



+ .3 % 



OO 



i d q 

X ad( ) Specific Volume Change 

-2X10"V° C -1.6X10" 3 

-8X10“ 4 /%o -2.4X10” 4 

-4 A 

-4X10 /lOOdb -4X10 



Pressure 



+100 db 



It has already been shown that Internal waves are oscil- 
lations of isodensity surfaces. Since temperature is the 
dominant parameter within the thermocline, the motions of an 
isodensity surface may be approximated by those of an iso- 
therm. This approximation is very convenient since 
temperature is much easier to measure than density. Temper- 
ature measurements for this reason have been widely used for 
obtaining observational data. 



APPENDIX B 



Interpolation Scheme 

The following parameters were sequentially sampled 
from each sea-sensing unit: 

D a * *^i » *^2 * t 3» T 4 * 
where D = depth and T = temperature. 

Then the depth of each thermometer when sampled 

may be found to first order by 

ADi = ( D b~ D a ) 1 * 1 = 

5 ( 12 ) 

D 1 = D a + A D i + Si 

S t h 
^ Is the vertical offset of the i thermometer from 

the servo-system thermistor. 

The depth of the Isotherm ( Drp ) to first order Is 

then found by using the two temperature measurements which 

bracket the Isotherm: 

D t = D a + [ (Di+l - ©i)/(Tl.--T 1+1 )j (T A - T) 

where T Is the temperature of the desired isotherm. 
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APPENDIX C 



Doppler Effect 

Phase velocity, c, Is defined as 
C = w/k = f/X' 1 



When sampling at three points simultaneously, the Doppler 
shift occurs only In frequency, leaving the wavelength 
unaffected. Thus: 



C +AC 



GP+ A CD 
k 



f +Af 
X' 1 



Ac = Af /X -1 

For frequency 9.47 cph of series 8, X ,== 8.62 X 10 3 ITI 1 . 
X _l = 8.62 X I0~ 3 m occurs at 7.90 cph on the first modal 
curve of Figure 6. 



Ac = 



9.47-7.90 
8.62 x I0" 3 



x iq2 

3600 



1.57 ^ e Cm 

, 3 i " sec 
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